function expresult = eigqr(A)
    [V,D]=eig(A)
    [Q,R]=qr(V)
    values = diag(D);
    exps = exp(values);
    Y = R \ Q'
    vl = V*diag(exps)
    % hint uit help en opgave : A/B = (B'\A')'
    expresult = (V*diag(exps))*Y;    
end